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a. 


Consider  Solutions  <G(x,e),H(x,E)> 


of  the  von  Karman  equations  for  the 


swirling  flow  between  two  rotating  coaxial  disks 


(1.1  ) 


EH 


IV 


+  HH' • •  +  GG* 


and 

(1.2) 

with  boundary  conditions 


EG"  +  HG"  -  H'G  =  0 


V,  '/  s'  \  ^  V 


(1.3)  H(0,£)  =  H' (0,€)  =  H( 1 ;e)  =  H' (1 ,£)  =  0 

(1.4)  ^  G(0,£)  =  s,  G(1,£)  =  1,  |sl  <  1  .  , 

In  this  work  we  establish  the  existence  of  solutions  for  t  small  enough.  In 
fact,  if  n  is  a  given  positive  integer  with  sign  s  =  (-1)"^  then  there  is  - 
for  '£  >  0  sufficiently  small  -  a  solution  with  the  additional  property: 
G(xr^)  has  n  interior  zeros.  If  n  >  1  there  are  at  least  two  such 
solutions.  If  s  =  0  there  is  at, least  one  such  solution  for  every  positive 
integer  n.  The  asymptotic  ’^shape'^^f  these  solutions  is  described.^ 


AMS  (MOS)  Subject  Classifications:  34B15,  34E15,  35Q10 
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SIGNIFICANCE  AND  EXPLANATION 


Under  appropriate  conditions  the  steady-state  flow  of  fluid  between  two 
planes  rotatinq  about  a  common  axis  perpendicular  to  them  may  be  described  by 
two  functions  H(x,e),  G(x,e)  which  satisfy  the  coupled  system  of  ordinary 
differential  equations 

+  HH* • •  +  GG*  =0 
CG"  +  HG*  -  H'G  =  0  . 

Ibe  quantity  e  >  0  is  related  to  the  kinematic  viscosity  and 
usually  called  the  Reynolds  number. 

ITiese  equations  have  received  quite  a  bit  of  attention.  First  of  all, 
people  who  are  truly  interested  in  the  phenomena  modeled  by  these  equations, 
e.q.  fluid  dynamicists,  are  interested  in  this  problem.  Howeyer,  as  these 
equations  have  been  studied  by  a  variety  of  mathematical  methods,  they  have 
taken  on  an  independent  interest.  The  major  methods  employed  have  been  (i) 
Matched  Asymptotic  Expansions  and  (ii)  Numerical  Computations.  In  both 
approaches  technical  problems  have  appeared.  There  may  be  "turninq  points," 
i.e.  points  at  which  H(x,€)  =  0.  Such  points  require  special  and  delicate 
analysis  within  the  theory  of  (i).  As  numerical  problems,  these  equations  are 
"stiff"  -  precisely  because  e  is  small.  The  occurrence  of  "turninq  points" 
only  makes  computation  more  difficul't. 

For  these  reasons,  these  equations  have  become  "test"  problems  for 
methods  of  "matchinq  in  the  presence  of  turninq  points"  and  "stiff  O.D.E. 
solvers."  However,  when  one  has  "test  problems,"  one  needs  to  know  the 
answers.  Unfortunately  here  the  answers  are  larqely  unknown. 

In  this  report  we  establish  the  existence  of  a  multiplicity  of  solutions 
of  small  e  >  0.  In  addition  we  qive  a  complete  asymptotic  description  of  the 
shape  and  structure  of  these  solutions.  These  results  confirm  the  formal 
asymptotic  work  of  Watts  [241  and  qo  much  further.  These  results  should  be 
useful  in  checkinq  numerical  studies  of  this  problem. 


The  responsibility  for  the  wordinq  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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Introduction 

In  1921  T.  von  Karman  [5]  developed  the  similarity  equations  for  axi-symmetric, 
incompressible,  steady  flow  -  "swirling  flow*.  Let  (q^,q0»qjj)  be  the  coordinates  of 
velocity  in  cylindrical  coordinates,  (r,9,x).  von  Kaman  assumed  that  there  is  a  function 
H(x,£)  such  that 

=  -H(x,e)  . 

Then,  as  a  direct  consequence  of  the  steady  state  Havier-Stokes  equations  one  finds  that 
(see  (11,  (5])  there  is  a  function  G(x,e)  so  that  the  velocity  components  are  described  by 

=  “I  H' (x,e),  qQ»'|G(x,e)  . 

The  functions  < G(x,e ),H(x,e )>  satisfy  the  equations 

(1.1)  +  HH' •  '  +  GG‘  =  0  , 

(1.2)  EG"  +  HG*  -  H’G  «  0  . 

The  quantity  e  >  o  is  related  to  the  bulk  viscosity.  Equation  (1.1)  can  be  integrated  to 
yield 

(1.3)  eH'”  ♦  HH*  +  G^  -  j  (H')^  »  U 
where  p  is  a  constant  of  integration. 

In  the  case  originally  studied  by  von  Karman,  the  flow  above  a  single  disk,  we  have  a 
problem  on  the  infinite  interval  [0,»1  and  the  constant  of  integration  is  known,  i.e., 

where  Sl^  =  G(“,€). 
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If  wa  conaldar  th«  flow  batwaan  two  planaa,  x  •  0i  x  •  1  rotating  with  constant 
angular  valocltlas  Qq/2«  0^/2,  than  tha  quantity  u  >  U(e)  la  unknown.  This  lattar  case 
was  first  studlad  by  Batchalor  tl]  and  Stawartson  [22]  who  gava  conflicting  arguaents  and 
conjecturaa.  In  this  casa  tha  boundary  conditions  ara 
(1.4b)  H(0>e)  «  H(1.e)  •  0,  [no  panatratlon] 

(1.4b)  H'(0,e)  •  H'd.e)  >  0,  [no  slip] 

(1.4e)  G(0,e)  -  0^,  Gd.e)  -  0^,  l%l  +10,1*0. 

Many  authors  hava  studlad  this  problaa.  Nuaarlcal  calculations  hava  baen  carrlad  out 
by  Lanca  and  Rogara  [9].  C.  E.  Paarson  [IS].  D.  Graanspan  [3],  D.  Schults  and  D.  Graenspan 
[19],  L.  0.  Wilson  and  N.  L.  Schryar  [26],  G.  L.  Nallor,  P.  J.  (Siappla  and  V.  K.  Stokes 
[12],  N.  D.  Nguyen,  J.  P.  Rlbault  and  P.  Florent  [13],  S.  M.  Roberts  and  J.  S.  Shipman 
[IS],  H.  B.  Ksllar  and  R.  K-H.  Szeto  [6].  Formal  matched  asymptotic  expansion  methods  have 
been  used  by  A.  Watts  [24]  (who  also  did  numerical  calculations),  K.  K.  Tam  [23], 

H.  Rasmussen  [17],  B.  J.  Matkowsky  and  W.  U  Slegmann  [11].  Undoubtably  many  others  hava 
also  worked  on  this  problem. 

Nevertheless,  the  basic  questions  of  "existence*  and  "uniqueness*  hava  ramalned 
unanswarad.  S.  P.  Hastings  [4]  and  A.  R.  Elcrat  [2]  have  proven  exlstance  and  uniqueness 
for  large  e,  nieir  arguments  are  essentially  a  perturbation  about  e  ■  *.  j.  B.  McLeod 
and  S.  V.  Barter  [10]  considered  the  special  case  where  *  0.  They  have  shown  the 

existence  of  a  solution  for  all  e  >  0  and*  for  those  solutions,  they  gave  a  complete 
discussion  of  the  asymptotic  behavior.  These  solutions  are  bounded.  In  fact 

H(x,c)  .  0(/e) 

G(x,c)  -  0(1)  . 

There  have  been  no  results  on  the  number  of  solutions.  In  [21]  J.  Serrin  referred  to 
computational  results  of  C.  E.  Pearson  in  tha  case  ■  -R^  *  0.  Pearson's  computations 
indicated  the  existence  of  solutions  which  are  not  anti -symmetric  about  x  -  -  -j.  Thus, 
since  <G(1  -  x,e),-H(1  -  x,e))  is  a  solution  whenever  (G(x,e),H(x,e) >  is  a  solution,  these 
computational  results  lead  one  to  conjecture  that  there  ara  several  solutions.  Indeed  much 
of  the  computational  work  since  Pearson's  early  paper  seems  to  confirm  this  conjecture. 
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Introduction 

In  1921  T.  von  Karman  [5]  developed  the  similarity  equations  for  axi-symmetric, 
incompressible,  steady  flow  -  "swirling  flow".  Let  (q^,qg,qjj)  be  the  coordinates  of 
velocity  in  cylindrical  coordinates,  (r,6,x).  von  Kaman  assumed  that  there  is  a  function 
H(x,E)  such  that 

=•  -H(x,e)  . 

Then,  as  a  direct  consequence  of  the  steady  state  Navier-Stokes  equations  one  finds  that 
(see  [1],  (51)  there  is  a  function  G(x,e)  so  that  the  velocity  components  are  described  by 

q^  “jH'Cx.e),  qg  »  i.  G(x,e)  . 

The  functions  (  G(x,e ),H(x,e ) )  satisfy  the  equations 


(1 

.1 ) 

+  HH'  • 

•  +  GG* 

(1 

.2) 

EG" 

+  HG’ 

-  H’G  =• 

The  quantity  e  >  o  is  related  to  the  built  viscosity.  Equation  (1.1)  can  be  integrated  to 
yield 

(1.3)  €H*"  +  HH*  +  -j  G^  -  (H’ )^  »  M 

where  p  is  a  constant  of  integration. 

In  the  case  originally  studied  by  von  Karm&n,  the  flow  above  a  single  disk,  we  have  a 

problem  on  the  infinite  interval  lO,")  and  the  constant  of  integration  is  known,  i.e., 

M  =  ^8^ 

2  » 

where  8^  =  G(",e). 
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If  wa  conaldar  th«  flow  batwaan  two  planaSf  x  •  0>  x  ■  1  rotating  with  conatant 
angular  valocltlaa  Ojj/2,  Q^/2,  than  tha  quantity  p  -  p(e)  ia  unicnown.  Thla  lattar  eaaa 
waa  firat  atudlad  by  Batchalor  [1]  and  Stawartaon  [22]  who  gava  conflicting  argunanta  and 
conjacturaa.  In  thia  caaa  tha  boundary  condltlona  are 
(l*4a)  H(0,e)  •  H(l,e)  •  o,  [no  panatration] 

<l.4b)  H>(0,e)  >  a‘(1,e)  >  0,  [no  alip] 

<i.4c)  aio,€)mU^,  G(i,c)  -  a^,  la^jl  +  |aj  a  0  . 

Many  authora  hava  atudlad  thia  problaa.  Nuaarical  calculatlona  hava  baan  carriad  out 
by  Lanca  and  Rogara  [9],  C.  B.  Paaraon  [IS],  D.  Graanapan  [3],  D.  Schultz  and  D.  Graanapan 
[19],  L.  O.  wilaon  and  N.  L.  Schryar  [26],  G.  L.  Hellor,  P.  J.  Chappla  and  V.  K.  Stokaa 
[12],  N.  D.  Nguyan,  J.  P.  Rlbault  and  P.  Florant  [13],  S.  M.  Roberta  and  J.  S.  Shlpaan 
[18],  H.  B.  Kallar  and  R.  K-H.  Szeto  [6].  Poratal  aatchad  aayaptotic  axpanaion  nathoda  hava 
baan  uaed  by  A.  Watta  [24]  (who  alao  did  nuaarical  calculatlona),  K.  K.  Tan  [23], 

H.  Raaauasan  [17],  B«  J,  Matkowaky  and  w,  i,.  Slagnann  [11].  Undoubtably  nany  othara  hava 
alao  worked  on  thia  problan. 

Navarthaleaa,  tha  baaic  quaationa  of  "axiatanea''  and  "uniquannaa”  hava  ranainad 
unanawarad.  S.  P.  Haatlnga  [4]  and  A.  R.  Blcrat  [2]  hava  proven  exiatanca  and  uniquanaaa 
for  large  e,  Thalr  argunanta  am  eaaentially  a  perturbation  about  e  >  •>.  J.  B.  McLeod 
and  S.  V.  Parter  [10]  conaidared  tha  apacial  caaa  where  >  -Q^  *  0.  They  hava  ahown  tha 
exiatanca  of  a  aolution  for  all  e  >  0  andi  for  theaa  aolutlona,  they  gave  a  conplete 
diacuaaion  of  tha  aayaptotic  behavior.  Iheaa  aolutlona  are  bounded.  In  fact 

H(x,e)  -  0(/F) 

G(x,e)  -  0(1)  . 

Thera  have  baan  no  reaulta  on  tha  number  of  aolutlona.  In  [21]  J.  Sarrin  referred  to 
conputational  reaulta  of  C.  E.  Pearaon  in  the  caaa  Rq  •  P  0.  Paaraon'a  conputationa 
indicated  the  exiatanca  of  aolutlona  which  are  not  antl-aynnatric  about  x  -  -  Thua, 
alnca  <G(1  -  x,e),-H(1  -  x,e)>  ia  a  aolution  whanavar  <G(x,e),H(x,e) )  ia  a  aolution,  theaa 
conputational  raaulta  lead  one  to  conjecture  that  there  are  aeveral  aolutlona.  Indeed  much 
of  tha  conputational  work  ainca  Paaraon'a  early  paper  aaena  to  confirn  thia  coi;jacture. 
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Furthermore,  the  physical  boundary  conditions  (1.4a),  (1.4b)  together  with  the  conditions 
(I.Sa),  (I.Sb)  imply  that  (in  the  limit  as  e  ^  0+)  the  function  H(x,e)  may  not  have 
nodal  zeros,  i.e. 

(1.6)  Um  H(x,e)  nay  not  change  sign  . 

Therefore  w  turned  to  the  question  of  the  existence  of  "pathological*  solutions 

<G,H)  which  satisfy  the  differential  equations  (1.1),  (1.2),  the  bounds  (I.Sa),  d.Sb), 

the  boundary  condition  (1.4a),  (1.4b)  and  also  satisfy  (1.4c)  with 

(t.7a)  |no<e)|  +  |0,<  =  >l  *  0>  0  <  e 

(1.7b)  Lim  {|a  (e)|  +  |n  (e)|}  -  o  . 

€♦0 

Moreover,  (1.6)  implies  these  solutions  should  be  essentially  "positive"  in  the  sense  that 

H(x)  •  Lim  H(x,e)  >  0  . 

RemarX;  The  hypothesis  (1.5a),  (I.Sb)  implies  that  -  after  the  extraction  of  a  subsequence 

if  necessary,  the  functions  H(x,e)  are  convergent  to  a  non-trivial  limit  function. 

The  results  of  (8)  imply  that  H(x)  must  then  have  a  special  form.  That  is,  there 

are  n  numbers,  o  , a.,..,, o  ,  with 
0  1  n 


(1.8a) 

And,  on  the  Interval 


0  -  Oq  <  0i  < 


<0  »  1 


j  *  Opifeeeyn  -  1 
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i 


we  have 


(1.11)  H'(Oj(e),e)  «  0,  H"(o^(e),e)  >  o  . 

Hor cover, 

H(o^(e),e)  > 

and,  between  the  o^(e),  H(x,e)  Is  essentially  positive.  In  fact,  H(x,e)  is  essentially 

given  by  {1.8b)  while  the  function  G(x,e)  has  the  fora 

(1.12a)  G(x,e)  -  (-1  )^T^H(x,e),  o^(e)  <  x  <  0j^,(6) 
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where 


(1,12b)  T.»  '  —  ■  ■  -  j  *  0,1f#,*#n  -  1  • 

1  o.  ,  -  o. 

1+1  i 

Finally,  there  is  a  constant  t  <  n  such  that 

(1.13)  Lira  T.(e)  =  |tP“^  . 

£+0  ^ 


Remark : 

The 

characterization  of 

G(x,e)  qiven  by  (1.2a),  (1.2b)  can  be 

made  more  precise. 

Case  1 : 

‘>0  ^ 

0  . 

In  this 

case 

G(x,E) 

has  exactly 

n  interior 

zeros,  say  Y,,Y-,...,Y  , 

1  2  n 

and 

(1.14a) 

I.  <  0., 

=  0(eV3,. 

( 1 . 1 4h ) 

-  Y.  * 

Case  2: 

1  1 

‘’o  < 

0  . 

In  this 

case 

G(x,e) 

has  exactly 

(n  +  1)  zeros,  say  Y, ,Y,, . . . ,Y 

u  1  n 

and  ( 1 . 1 4a  ) , 

(1.14b) 

hold. 

Moreover 

(1.14c) 

‘’o  ' 

°  <  Y„ 

-  Op  .  0(e’^^)  . 

Once  one  has  proven  this  result  we  obtain  the  patholoqical  solutions  on  the  interval 
(0,11  by  takinq 


(1.15a)  Z  =  x/l,(e),  £•  =  e/L(e)  , 

(1.15b)  H(x,R’)  =  H(x,e),  G(x,e’)  =  l(e)G(x,e)  . 

The  functions  <G,H>  satisfy  (1.1),  (1.2)  -with  x  replaced  by  x  and  E  replaced  by 
E'.  In  addition  these  functions  satisfy  the  boundary  conditions  (1.4a),  (1.4h),  (1.4c) 


with 


Finally,  if  we  set 
( 1 . 16a  ) 


"o  =  V®’’  =  , 

n,  =  «,(£•)  =  q^L^^^(EM^^^  . 

e* 

^  =-rTrTe-^=  ■ 


(1.16c)  G(x,e)  -  G<x,e*)  -  0(e"^) 

we  obtain  a  solution  of  (1.1),  (1.2),  (1.4a),  (1.4b),  (1.4c)  with  e  replaced  by  e,  x 
replaced  by  x,  with  (1^  replaced  by  0^,  and  Sl^  replaced  by  £1^  where 

“o”V^1'  =  1  • 

For  completeness  sa)ce  we  formulate  this  last  result  as 
Theorem  II;  Let  n  >  1  be  an  integer.  Let  s  be  a  given  real  number.  Then  for  e 
small  enough  there  is  a  solution  (  G,H  >  of  (1.1),  (1.2),  (1.4a),  (1.4i  ),  (1,4c)  with 
(1.17)  “o  “  =  1  . 

This  solution  may  be  described  in  a  manner  similar  to  the  description  given  in  Theorem  1. 
There  are  exactly  (n  +  1)  numbers 

0  =  <  0,(E)  <  ...  <  0^_,(E)  <  0^(6)  =  1 

at  which  H(x,c)  has  its  relative  minima,  i.e.. 


(1.18a) 


H'(Oj(e),e)  =  0,  H"(a^(e),e)  >  o  . 


Moreover,  between  the  o^(e)  the  function  H(x,E)  is  essentially  positive.  That  is,  for 
any  given  6  >  0,  26  <  we  have,  for  small  e. 


(l.iab) 


Furthermore 


H(x,e)  >  0,  0^(6)  +  6  <  X  <  -  6  . 


(E)-'!  yc))  =  S  = 


The  function  G(x,£)  has  at  least  n  nodal  zeros;  0  <  Y.tE)  t  Y  (e)  <  ...  <  y  (e)  <  1. 

1  2  n 

Moreover 

(1.20a)  Yj(e)  <  o^(£),  '  0(E),  j  =  1,2,...,n  . 

If  s  ^  0  and 

sign  8  =  (-1)1''*'^ 

then  G(x,E)  has  (n  +  1 )  zeros.  The  additional  zero,  7p(E)  satisfies 
(1.20b)  0  <  Y^j(E)  =  0(e)  . 

Furthermore 

(1-21  )  BHl  «  Bon  ^  (e"^)  , 
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Renarit:  One  can  choose  to  characterise  the  solution  (  G,H  >  by  the  number  of  "humps*  or  by 
the  number  of  interior  (nodal)  zeros  of  G(x«e)>  Suppose  we  choose  to  discuss  the  number 
of  humps.  If  s  ^  0  let  <G(x,e).H(x.c))  be  the  solution  described  In  Theorem  II  with 
n  humps  and 

i  -  1,  e  -  |a|e  . 

Then 

H<x,e)  .  -  H(1  -  x,e)  , 

G(x,e)  -  j  G(1  -  x,i) 

Is  another  solution  of  (1.1),  (1.2),  (1.4a),  (1.4b)  and  (1.17).  On  the  other  hand.  If  one 
chooses  to  look  at  the  number  of  Interior  zeros  of  G(x,e)  we  have  the  following  situation 
Case  1 ;  s  >  0.  For  every  even  n  >  2  there  are  at  least  two  solutions  (G,H>,  (G,B>  of 
(1.1),  (1.2),  (1.4a),  (1.4b)  and  (1.17)  with  6(x,e),G(x,e)  having  exactly  n  Interior 

zeros  and  which  also  satisfy 

H(x,e)  >  0,  H(x,e)  >  0,  (essentially) 

Let  nan  and  <G(x,e),H(x,e)>  be  the  solution  described  In  Ttieorem  II.  From  (1.9)  we 
see  that 

g,  >  0,  go  >  0  • 

Hence  G(x,e)  has  exactly  n  «  n  Interior  zeros.  In  addition,  let  n  •  n  -  1  and 
<G(x,e),H(x,e)>  be  the  solution  described  In  Theorem  II.  Then 

g^  <  0,  g^  <  0  (essentially) 
and  G(x,c)  has  exactly  n  t  1  -  n  Interior  zeros. 

Case  2i  s  <  0.  For  every  odd  n  >  3  there  are  at  least  two  solutions  <G,H>,(G,H>  of 
(1.1),  (1.2),  (1.4b),  (1.4b)  and  (1.17)  with  G(x,e),G(x,e)  having  exactly  n  Interior 

zeros  and  also  satisfy 

M(x,e)  >  0,  H(x,e)  >  0,  (essentially)  . 

If  n  a  1  there  Is  at  least  one  solution  <G,H>  of  (1.1),  (1.2),  (1.4c),  (1.4b),  and 
(1.7)  with  G(x,e)  having  exactly  n  a  i  Interior  zeros  while 

H(x,e)  >  0  (essentially). 

Let  n  a  n  and  ( G(x,e),H(x,e))  be  the  solution  described  In  nieorem  II.  From  (1.9)  we 
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see  thdt 


g,  <  0,  go  >  0  • 

Hence  U(x,e)  has  exactly  n  n  interior  zeros.  If  n  >  1  let  n  »  n  -  1  and  let 
<G(x,e),H(x,e))  be  the  solution  described  in  Theorem  II.  Then 

q,  >  0,  go  <  0 

and  G(X(C) 

Case  3;  8=0.  For  every  n  >  1  (even  or  odd)  tnere  is  at  least  one  solution 

<G(x,e),H(x,e))  with  G(x,e)  having  exactly  n  interior  zeros  and 

H(x,e)  >  0  (essentially)  . 

Let  n  =  n  and  let  ( G(x,e ) ,H(x,e)>  be  the  solution  described  in  Theorem  II. 

The  basic  Theorem  I  is  proven  via  a  ‘'shooting*  argument.  The  basic  estimates  follow 
from  the  following  analysis.  When  H(x.c)  is  small,  i.e. 

H(x,e)  =  0(E^/^) 

then  one  studies  the  "stretched"  problems  let 


(1.22a) 


sV3  ' 


h(C,E)  =  e'^/^H(x,E),  g(5,E)  =  e"^/^G(x,E)  . 


h"‘  +  hh"  +  J  E^/^g^  _  1  (h')^  =  =  Ue’^^ 

g“  +  hg'  -  )>’g  =  0  . 


(1.22b) 

The  functions  <g,h>  satisfy  the  equations 
(1.23a)  h"  ■  +  hh"  +  —  £-•  g 

(1.23b) 

With  M  =  0(e),  i.e.,  u  =  0(1).  We  find  that 

(1.24a)  h(C.E)  ♦  h(C)  ;  a  quadratic  function 

of  the  form 
(1.24b) 

Furthermore 

g(t.e)  ♦  g(C) 

where  g(C)  satisfies 

(1.25)  g"+ng'-R'g=0. 


h(e)  =  I  (C  -  Cq)^  . 
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The  eolutlona  of  thla  problem  are  discussed  in  the  Appendix.  On  the  other  hand,  when 
H(x,e)  is  “lar^e*.  then  the  development  in  18)  shows  that  H,H' ,H* ,H' ' ' ,G,G' .ti*  can  all 
be  estimated  in  teres  of 

|G<x.e)| 

|H(x,e)|  • 

Fortunately,  we  do  not  require  that  H(x,e)  be  too  large.  In  fact, 

H(x,e)  > 

is  sufficient.  Hence  the  requirements  of  “small*  H{x,£)  and  “large*  H(x, e)  overlap  and 
we  are  able  to  give  a  complete  analysis. 

Realising  these  facts  one  proceeds  as  follows. 

Starting  Procedure  (See  Theorem  3.2.) 

For  every  choice  of  M.h^.a.g^  there  is  a  solution  (g.h)  of  (1.23a),  (1.23b)  on  the 
interval  (0,a]  which  also  satisfies  the  boundary  conditions 
(1.26a)  h(0,e)  -  h'(0,e)  •  0,  h“(0,t)  -  hj  >  0  , 

(1.26b)  g(0,e)  -  g^,  g(a,e)  ■  h(o,e)  . 

The  results  of  Section  3  show  that  this  solution  (G,h>  of  (1.1),  (1.2)  -  originally 
defined  only  on  the  Interval  may  be  continued  to  entire  interval  (0,2s  -  5) 

and,  on  this  interval 

“2 

(1.27)  H(x,C)  ■  —  (1  -  cos  Tx),  G(X,e)  »  TH(X,E) 

where  0  <  t  <  1  and  t  1  as  a  *  The  results  of  Section  4  show  that  this  solution 

may  in  fact  be  continued  to  the  larger  interval  2s[i  +  — ^ — ]  -  6  where  t  is  a  negative 

|7| 

number  described  in  the  Appendix.  Furthermore,  for  small  e  we  have 

h 

{1.28a)  H(x,e)  "  (1  -  cos  t  (x  -  2s)),  2s  +  6  <  x  <  2s[l  +  — ]  -  5  , 

X?  1^1 

(1,2bb)  G(x,e)  •  -T,H(x,c),  2s  +  6  4  X  <  2sli  +  —2—]  -  4  , 

(1.28c)  •  t|t|  . 

Thus  we  have  exhibited  2  "humps*.  Proceeding  in  this  way  we  construct  a  solution  with 
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n  hunps.  liet  - 
employ  an  elementary 

provided  that 


X  (e)  be  the  n'th  relative  minimum  of  H(x>e).  In  Section  5  we 
n 

degree  theory  argument  to  show  that  one  may  choose  Urh^  so  that 

H(x^,e)  -  0  , 
n 

G(x^,e)  > 

sign  q,  -  (-I)"*’  . 


In  this  way  we  prove  Theorem  I< 


2,  Existence  of  Solutions  Away  from  Turning  Points 
In  this  section  we  are  concerned  with  solutions  of  the  equations  (1.1),  (1.2) 

(2.1)  +  HH"'  +  GG’  -  0  , 

(2.2)  eG“  +  HG*  -  H'G  «  0  , 
with  initial  data 

(2.3a)  d'’H(XQ)/dx'’  »  v  »  0,1, 2, 3,  d'’G(Xg)/dx''  «  G^,  v  «  0,1  , 

where 

(2.3b)  Hjj  >  0  . 

In  dealing  with  this  case  we  use  a  basic  estimate  of  Kreiss  (lemma  2.1  of  [7])  which  we 
include  for  the  sa)ce  of  completeness. 

Lemma  2.1;  Consider  the  differential  equation 

(2.4)  ®  5x  “  F(x),  o  <  X  <  8  , 

where  a,  F  are  continuous  functions  with 

a(x)  >  0 

and  e  >  0  is  a  positive  constant.  The  solutions  of  (2.4)  satisfy  the  estimate 

(2.5)  |y(x)|  <  iF/al  +  0(x,o) |y(a) I ,  a  <  x  . 

Here 

(2.6a)  ®  U(t)( 

'*  a<t<x 

(2.6b)  o(x,(>)  »  exp{-  —  /  a(t)dt)  . 

a 

Proof !  The  solutions  of  (2.4)  are  given  by 


y(x)  »  y(a)o(x,a) 


+  /  ex  Aj  /  a(s)dB} 
a  X 


We  rewrite  this  as 


y(x) 


y(a)o(x,a) 


a(t ) 
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For  fixed  x  >  a  the  function 


a(s  )ds) 


is  monotone  increasing  as  t  increases  from  a  to  x. 
the  mean  value  theorem. 

Lemma  2.2;  Let  <  G(x,e  ) , H(x,e  )>  be  a  solution  of  (2.1) 


The  estimate  (2.5)  follows  from 


,  (2.2),  (2.3a),  (2.3b). 


Let 


1 


1  * 

s(x,Xq)  =  expl-  —  /  H(t,e)dt) 


Let 


G"(Xq,£)  =  Gj  . 


Then,  for  x  >  Xy  we  have 


<2. 7a) 


|G'(x,«)|  <  IG/HI  •  IH'I  +  s(x,Xn)|G, I  , 

*0'  0' 


(2.7b)  |G"(x,e)l  <  IG/HI^  ^  •  IH"I  +  S(X,X  )|G J  , 

X^^X  U  A 


(2.7c)  |H***(x.e)|  <  IG/Hl  •  IG*I  ♦  s(x,Xrt))H,)  . 


Proof ;  The  estimates  (2.7a),  (2.7c)  follow  from  Lemma  2.1  and  equations  (2.2),  (2.1) 

respectively.  Differentiating  (2.2)  we  have 

(2.b)  ec' • •  +  HG"  »  H"G  . 

The  estimate  (2.7b)  follows  from  (2.B)  and  Lemma  2.1. 

Lemma  2.3:  Let  G,H  be  a  solution  of  the  above  problem  in  some  interval  Xq  <  x  <  x^ 
with  the  following  properties 

(2.9)  H  >  0,  IG/Hl  <  H,  X,  -  Xq  <  min(1/M,1)  . 

*0'*1 


Then  there  are  constants  which  deijend  only  on  ri^,  V  =  0,1, 2, 3;  G^,  v  =  0,1  and 

M  and  nut  on  c  such  that 
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Id-'G/dx-’l 


j  -  0,1 


(2.10) 


<  K 


Ij' 


id^H/dx^i 


3  »  0,1 ,2,3  . 


Also 

(2.11)  IG'I^  *'"‘22  *  l°2*'  °2  "  “"<*0'^’  • 

Proof!  For  any  solution  of  the  above  equations  we  may  apply  Letana  2.2  and  obtain  the 
estimates  {2.7a),  {2.7b),  (2.7c).  Also,  the  Taylor  expansion 


H’(x)  -  ♦  (X  -  Xq)H2  +  /  (6  -  XQ)H”'(C)dC 


qives  us 
(2.12) 


IH'I  <  |H,|  ♦  |x  -  XqIHj  +-1  (X  -  *o>^'“’"'x„,, 

0  y 


Therefore  by  (2.9)  and  (2.7) 

.2 


IH"  •  I  <  h' 

“o'*! 


IH’I  +  ►‘|G,1  ♦  IH3I  <  J  (X  -  X  )VlH-"l  + 

*0'  1  J  X  1  U  Xq,X^ 


+  +  (X,  -  XjjjlHjj}  +  m|gJ  +  |Hj1 


By  (2.9)  (x,  -  Xq)V  <  1  and  x,  -  Xq  <  1.  Therefore 

lH'”l  <  2m2(|h,|  +  IHjI)  +  2h|gJ  +  2IH3I 


Thus  we  have  proved  the  estimate  for  H'".  By  Taylor  expansion  (see  (2.12))  we  obtain  the 
estimates  for  and  by  (2.7a)  and  (2.7b)  they  follow  also  for  G,G',G". 

We  shall  now  use  these  estimates  to  derive  existence  theorems. 

Theorem  2.1:  Consider  the  initial  value  problem  (2.1)-(2.3)  and  assume  that 

(2.13)  Hjj  >  6  >  0  . 

Then  there  is  an  Interval  Xq  <  x  <  Xj,  x^  -  Xq  >  0,  independent  of  e,  in  which  the 
above  problem  has  for  all  e  with  0<e<1  a  unique  solution.  Moreover  the  estimates 


(2.10)  of  Lemma  2.3  are  valid.  The  constants  depend  only  on  v  >  0, 1,2,3; 

G^,v  »  0,1,2  and  5. 

Proof ;  Let  e  >  O  be  fixed.  From  the  general  existence  theory  for  ordinary  differential 
equations  it  follows  that  there  is  an  interval  Xq  <  x  <  x^  where  the  conditions  of  Lemma 
2.3  are  satisfied  with  M  -  2(|Gq/Hq|  +1).  We  want  to  estimate  x^.  Taylor  expansion 
gives  us 

|h(x)  -  HqI  <  <x  -  Xq)K2,,  |g(x)  -  GqI  <  (x  -  Xo)K,,  . 

Therefore  h(x)  >  Hg,  |g(x)|  <  |Gg|  ♦  |Hq(  and  |h(x)/G(x)|  <  M  for 

0  <  X  -  Xg  <  min(1/K,,,Hg/(21C2i))  *  "  Xq  * 

Thus  the  solution  exists  in  this  interval  ^uul  the  theorem  is  proven. 

Now  consider  the  limit  process  e  ♦  0.  We  want  to  prove 
Theorem  2.2t  Assume  that  G..,H.  are  functions  of  e  with 

(2.14a)  lie  Hy  »  H  ,v  »  0,1,2,  Urn  Gg  «  Gg,  Hg  >  S  >  0  . 

e*0  e*0 


Assume  also  that 

(2.14b)  H2,G^,G2  are  uniformly  tx>unded  . 

Let  Xg,Xf  be  as  Theorem  2.1.  Then  the  solutions  of  the  Initial  value  problem  (2.1) -(2. 3) 
converge  on  any  interval  Xg+4<x<x^,  5>0  to  the  solution  of  the  reduced  problem 
(2.15a)  HH'"  +  GG*  -  0,  d'H/dx''|  “  V  =  0,1,2  , 

X-Xq 

(2.15b)  BS*  -  B’C  «=  0,  C(Xq)  =  Gg  . 

Proof ;  (2.14)  and  (2.10)  show  that 

EG"  «  0(e)  ♦  0  . 

Differentiating  (2.1)  we  obtain  for  y  =  H*'''  the  equation 

ey*  +  Hy  +  H'H'"  +  (G*)^  +  GG"  »  0  . 

Therefore  by  Lemma  2.1  and  (2.14) 


(2.16)  |y(x)|  <  ^ +  |y(Xp)  |s(x,Xo)  . 

Xg.x 

By  (2.14)  eH'*'(Xy)  is  bounded.  Therefore  y(Xg)R(x,Xg)  Is  hounded  for  x  >  Xg  +  6. 
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Theretore  ey  *  0  and  the  theorem  tollovrs  by  standard  compactness  arguments.  Tnis  proves 
the  theorem. 

It  is  easy  to  see  that  the  solution  oi  (2. IS)  has  the  form 


H(x)  “  Mq  +  -^  sinT(x  -  x^)  +  — I  (1  -  C08T(X  -  Xjj)) 

(2.17) 

G(x)  »  TH(x),  t  =  Gq/Hq  . 

Up  till  now  we  have  only  proved  the  convergence  to  G,H  in  the  interval 
Xq  <  X  «  Xj.  However,  we  obtain  immediately  uniform  convergence  in  any  interval 
Xq  <  X  <  x^  where  H(x)  >  a,  a  any  constant  >  0.  This  we  can  also  express  in  another 
way.  Let  x  <  Xq  <  x  be  the  first  points  to  Che  left  and  right  of  Xq  with 
H(x)  °  H(x)  o  0.  Then  we  can  prove  existence  of  solutions  <G,H>  of  (2.1),  (2,2)  and 
uniform  convergence  Co  (g,h)  in  any  interval  x  <  Xq  <  x  <  x^  <  x.  (Of  course,  if  we 
move  Xq  Chen  we  have  to  change  Che  initial  conditions  to  obtain  the  same  (G,H)). 

In  (81  we  proved  that  a  necessary  condition  to  obtain  an  order  one  solution  of  the 
rotating  disc  problem  is  that  H' (x)  •  H’(x)  «  0  i.e.  we  can  write  (G,H)  in  the  form 


(2.18) 


cosT(x  -  x) ),  G  =  tii. 


X  <  X  <  X 


Hence,  we  shall  seek  such  solutions 


For  later  purposes  %<e  write  the  equation  in  another  forin«  We  inteqrate  (2>1) 


and  get 

(2.19)  eH"'  ♦  HH"  ♦ (G^  -  h'2)  -  ti  . 

To  obtain  a  limit  solution  of  the  form  (2.18)  it  is  necessary  and  sufficient  that 
lim  u  *  0  because  a  simple  calculation  shows  that  (S^R)  satisfies  the  equation 

Bfl“  +  -j  (B^  -  (B*  )^)  »  0  . 

For  our  purposes  the  right  choice  is 

(2.20)  M  -  eti  . 

Instead  of  we  can  give  u  as  initial  condition  and  compute  from  (2.19). 

particular  is  bounded  if 

(2.21)  HyHj  +  (Gq  "  **1*  “  ^  bounded  . 

Finally  we  collect  a  number  of  formulas  which  we  will  need  later. 

Lemma  2.4i  We  can  write  the  equations  (2.1 )>  (2.2)  in  the  form 


(2.22) 


H"'(x)  »  H"'(Xg)e 


-i/\(n)dn  .l/\(n)dn 

*0  .  i  -  5 


G(5)G’(t)dC  , 


(2.23) 


GMx)  -  G'(XQ)e 


.  1  r  5 

+  -  /  e 

’‘o 


HMe)G(5)de  , 


(2.24) 

G'(x)  -  (G(x)/H(x))H’(x) 

-  eG"(x)/H(x) 

(2.25) 

^  (G/H)  =  (HG*  -  GH')/H^  • 

-eG"(x)/H^(x) 

Also 

(2.2b) 


G"(x) 


2 

e 


G" (XQ)e 


X 

I 


H(n)dn 


X 

/  H(n)dn 

^  H"(e)b{C)de  . 


Proof;  Those  equations  follow  directly  from  (2.1).  (2.2)  and  (2.8). 


In 
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Exiitance  of  Solutions  when  H*  is  Small 


We  conaider  again  the  initial  value  problem  (2.1)-(2.3)  and  write  (2.2)  in  the 
integrated  form  (2.13)  with  p  of  the  form  (2.14).  Consider  initial  data  of  the  form 

(3.1)  Hq  -  Hj  »  hj,  P  -  ep.  h^  >  0,  v-  0,1,2  , 

(3.2)  Gq  -  G,  »  9^  >  0,  V  -  0,1  . 

P,h^,g^  can  be  functions  of  e  but  we  assume  that 


U.g^j.g,  ,hQ,h^  ,h2  and 


h  h  I 

0  2  2  1  ^  r  .  i  gV3  2 

^1/3  *^2  ^0 


are  bounded  independently  of  e.  we  assume  also  that 

(3.3)  Gj  “  h,gQ  -  hgg,  >  0  , 
and  is  bounded  independently  of  c.  we  want  to  prove 

Theorem  3. 1 <  There  is  an  Interval  Xq  <  x  <  ,  x^  -  Xq  >  0  independent  of  e,  in  which 

the  above  problem  has  a  unique  solution  for  all  e  with  0  <  e  <  1 .  Moreover 

(3.4)  0  <  -J  Hj  <  H*(x)  <  y  Hj,  x^  <  X  <  x^  , 

and  the  estimates  (2.5)  of  I.eama  2.1  are  valid,  x^  -  Xq  and  the  constants  depend 

only  on  ■  O,1,2,3>g0/h0  and  G^,v  »  0,1,2.  Also 

(3.5)  eG"(x)  »  0(e),  h(x>  >  H2(x  -  X0)^  . 

Proof I  Let  X0  <  X  <  x^  be  the  largest  interval  satisfying  (3.4).  Then  (2.20),  (3.2), 

(3.3)  and  (3.4)  imply  that  G''(x)  >  0.  Therefore  by  (2.19) 

IG/Hl  <  G0/H0  -  g0/h0  -  M  . 

0'*1 

Thus  Lemma  2.1  sho%ra  that  d^H/dx'',  v  «  0,1,2,3|  d^G/dx'',  v  -  0,1,2  are  bounded  if  (3.4) 
holds  and  x^  -  X0  <  mln(1/M,1).  By  Taylor  expansion 

|h"(x)  -  H"(Xo)l  <  (X  -  x„)lH"’l 

o  X,Xy 


Therefore,  we  can  find  x^,  independent  of  e,  such  that  (3.4)  is  valid.  Then  (3.5) 
follows  fro,«  Lemma  2.1  and  the  assumption  that  h0  >  0,  h^  >  0  and  that  G2  is  hounded 
independently  of  e.  This  proves  the  theorem. 
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Now  we  consider  the  limiting  process  e  >  0.  Assume  that 

*  o»  ♦  0,  h^  ♦  h^  >  0,  ♦  0,  ♦  o,  gg/h^j  ■*■  >  o  . 

-  1/3  2 

^O'^O'^I'^I  ^  large.  However,  we  assume  that  |i,c  '  g^  are  uniformly  bounded 

that 

(3.6)  g,  =  (go/hy)h,  -  t/hg,  hf  »  2hoh2  +  2e’''^p 

where  t  >  0,p  are  also  uniformly  bounded.  Tnis  assumption  guarantees  that  we  have  bounds 
for 

Gj  =  t  >  0  and  Hj  =  p  +  M  -  -i-  e^^^g^  . 

By  Section  2  and  Theorem  3.1  <  G,H  )  converge  in  a  neighbourhood  of  to  a  solution 

<S,H)  of  the  reduced  equation.  <  S,H  >  is  of  the  form  (2.12)  with  x  Xq.  we  want  to 

derive  a  relation  between  and  T.  By  (2.19)  and  (3.5) 


G(x^)  G(Xjj) 
H(x^  )  *  **(*0^ 


G"(x) 

H^(X) 


dx 


with 


dx 


/•  ’  G"  f* 

€  J  —  dx  <  const,  c  J  - ^ 

Xq  H  Xq  (Hp  +  -J  H2(x 


\2,2 
Xq)  ) 


const.  2  ^ 

it 


dx 


4V2 


2  2 
) 


<  const. 


_e 

,3/2 


const. 

k3/2 


(3.7) 


G(X^  ) 
H(X^  ) 


G(Xq) 

«<V 


Now  let  e  ♦  0  then  G(x^)/H(x.|)  ♦  T  and  therefore 

(3.8)  T  »  gy/hg  +  Od/h^/'^)  -  T,  +  0(1/h^/2)  . 

We  consider  now  a  two-point  boundary  value  problem  for  the  equations  (2.1),  (2.13)  in 


an  interval 


0  <  X  <  Xq  “ 


a  «  const.  >  0  • 

'  "  '  "0  ' 

Tne  bounUaiy  conditionw  are 

„.o,.KMO,.o.  OU,>.V^V^>°- 

Me  want  to  prove 

/->  51  (2.13)  (3.9)  has  «  solution  witfi  the 

Theorem  3.2:  The  .boundary  value  probleo  ( 2. 2 ) .  ( ^.  1  3  )  -  (  3  > 

r  r*  G"  H  K'.H'.H'"  can  be  estimated  by  hj.  t, . 
tollowinq  properties.  G.G'.G  ,«  ,« 

Also,  G.(x,)  >  0.  G-(x„)  >  0  provided  a  is  bid  enouqh. 

Proof ;  Introduce  new  variables  by 


V3'  o  .  H  »  e'^''^h  . 

c  X*  G  »  c  9/  ” 


Then  the  above  equations  become 
(3.10) 


...  **  1  # 

h+lih+-^l^  0 


(3.11 ) 


5  hq  -  gh  “  0  , 

h(0)  .  h(0)  -  0.  “  ^0>  * 


we  can  solve  the  reduced  equations 

•  i,- 

(3.12) 


h’  +  hi:  -  ^  -  0'  •'(0)  =  h(0)  -  0.  hto)  -  ^2 


J3.,3,  5  1  R|  -  qS  -  0.  g(0J  -  90' 

The  solution  of  (3.10)  is 
(3.14) 


R  -  -j  h^x^  . 


introducing  this  expression  into  <3..3)  and  the  boundary  conditions  gives  «s 

vg-.o.  gCO).go.  q'(a)  -  • 

«y  the  Appendix  the  general  solution  ot  the  differential  equation  (3.15)  is  of 

fotm  ,,  ,  1/1^. 

. ,  ( 0 )  ♦  0  and 

where  ^,(t)  oecays  exponentially  tor  t  .  , 


^2'"’  1  lim^ 

lin.^-rT-=1.  lim— 5^*1.  2 

t-H.  t  X— 


1  . 


Theretore,  if  »  is  sutficiently  large,  then 
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i 


(3.17) 


‘'o  - 


X  '  1  ,  X  ' 
2  1  1 


■(>,  (0) 


(3.18)  g  >  U,  q  >  0,  g/fi  '  for  all  sufficiently  large  x  . 

By  a  standard  perturbation  argument  it  follows  tnat  tne  full  system  (3.8),  (3.9)  has  a 
solution  witn 

\i  ^vj  \j_  '«.\i  V  V—  '“'V.  1/3 

(3.19)  Bd  g/dx  -  d  g/dx  (^  ■*■  **^  h/dx  -  d  h/dx  ^  ^  “  0,^,2  • 

hIso  h  =0  implies 
( J.2U) 

We  return  now  to  the  original  variables*  By  (3*17) 


*  *  *  1/3 

I  n  11^  <  C,e 

0,a  3 


2/3  1/3  ***1/3 

H(x)  =  X  +  0(e),  11*  (x)  =  h^x  +  0(£  '  H»(x)  =  +  0(e  '  ),  H*  "(x)  =  h  /z  '  . 

(3*21  ) 

G  =  +  0(e),  G*  =  e’^'^q  +  o(e^^^),  G"  =  g  +  0(c’''^)  . 

This  proves  the  theorem. 

Assume  now  that  a,  h-^i  qy,  and  (i  are  fixed  and  let  e  ♦  0.  At  Xy  “  the 

conditions  of  Theorem  3.1  are  satisfied.  Thus  <G,H>  can  be  continued  and  converges  to  a 
solution  <S,H>  of  the  reduced  equation  (2.12).  Here 

(3.22)  h2  =  H"(0)  = 

1  2 

and  satisfy  the  relation  (3.8)  with  *^0  ”  '2  *  0(e)  i.e. 

(3.23)  T  =  +  0(1/0^)  +  0(e)  . 

We  summarize  the  result  in 

Theorem  3.3:  Consider  the  two-point  boundary  value  problem  (2*1),  (2*13),  (3*9)*  Assume 
that  h2i  90/  and  M  are  fixed  and  a  sufficiently  large  let  e  0.  Then  (G,H> 

converges  uniformly  in  any  interval  0<x<x^<x  to  a  solution  <S,R)  of  the  reduced 
equation  (2.12)  with  T  and  satisfyinq  the  relations  (3.22)  and  (3.23)  respectively. 
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4.  Existence  of  Solutions  Through  a  Turning  Point 
In  the  previous  sections  we  have  shoinm  that  we  can  construct  a  solution  of  (2.1), 

(2.2)  for  0ax4x<x^  <x«  which  for  e  +  0  converges  to  a  solution  of  the 

reduced  equations  (2.12).  For  simplicity  we  assume  that  T  >  0  otherwise  we  change  G 
to  -G.  Also,  Xj  can  be  arbitrarily  close  to  x.  We  shall  now  show  that  the  solutions 
of  (2.1),  (2.2)  can  be  continued  through  the  turning  point  x.  I,et  x^  <  x  be  a  point 
near  x  where 

H(x,)  >  0,  H'(X, )  <  0,  H'(x,)  '  >  0,  G(x,)  >  0,  G’(x, )  <  0,  G“ (x, )  >  0 

for  all  sufficiently  small  e.  We  want  to  show 

Lemma  4.1.  Consider  (2.1)  in  the  form  (2.13)  and  assume  that  |i  is  given  by  (2.14).  For 

sufficiently  small  e  there  is  a  point  x^  >  x^  with 

H(x,)  -  >  0,  H'(x,)  -  <  0,  4  <  H"(x,)  •  h,  <4  3,  , 

(4.1)  1/^  2222 

G(X2)  -  ^  >  0.  GMXj)  «  e  g,  <  0  • 

Here  hg  >  0  is  a  sufficiently  large  constant  and  h^.gg.g^  are  constants  which  depend 
only  on  hg  and  not  on  e.  Also 

(4.2)  H(x)  >  0,  H*(x)  <  0,  G(x)  >  0,  G‘(X)  <  0,  G"(x)  >  0 

and  H'''(x)  is  uniformly  bounded  and  can  be  estimated  independently  of  ti  in  the  whole 

interval  x^  <  x  <  X2. 

Proof >  There  are  two  possibilities. 

2/3 

1 )  H(x)  >  e  hg  for  all  x  >  x^.  We  want  to  show  that  there  must  be  a  point  Xj 
where  H(x)  has  a  minimum.  Let  x^  <  x  <  x^  be  an  interval  where  G(x)  >  0,  H"(x)  >  0. 
Then  by  (2.20)  also  G''(x)  >  0  and  (2.19)  gives  us 


(4.3) 


0  < 


G(x) 

H(x) 


G(x,  ) 

hoTJT 


G(x  ) 

‘  hTx^' 


<  x  <  . 

4 


Thus  by  Lemma  2.1  H' ' '  (x),G'' (x)  are  uniformly  bounded.  Therefore,  choosing  x, 
sufficiently  close  to  x  and  x^  -  x^  sufficiently  small  guarantees 
(4.4)  T  ^2  ^  H"(x)  <  -j  x^  <  X  <  x^  . 

Also,  in  the  same  way  as  in  bection  3 
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G** 

—  dx  <  const. 


e  /  - -^^2 — - -  <  const, 

X,  (H(x^)  ♦•2  H2<*  ■  N” 


l.e.,  for  sufficiently  large 


(4.5) 


1  fllll  <  £<xi  <  fill! 

2  H(x^ )  H{x)  H(x, )' 


i.e.  G(x)  >  0  . 


Therefore  we  can  find  always  an  interval  x^  4  x  <  with  the  above  properties  whose 
length  x^  -  x^  does  not  depend  on  e  and  Xj.  Choosing  x^  sufficiently  near  to  x 
makes  H'(Xf)  as  small  as  we  like  because  H.H*  converge  to  R.R'.  Therefore  (4.4) 
implies  that  there  must  be  a  point  Xj  e  (x^.x^)  with  hMxj)  »  0.  At  this  minimum  (2.13) 
gives  us 

(4.6)  HH"  +  -i  »  eil  -  EH'  •  •  . 

For  sufficiently  small  e  and  x  -  x^  (4.5)  and  (2.6c)  show  that  H'*'  is  bounded. 

Therefore  H(x2)  *■  0(e)  which  is  a  contradiction. 

2/3 

2).  There  is  a  point  X2  with  h(x2)  “  ®  and  H'(x)  <  0  for  x,  <  x  <  X2. 

Using  (4.3)  we  find  again  that  G"(x)  >  0  and  that  (4.4)  and  (4.5)  hold  for 

2/3 

4  X  <  X2>  In  particular  G(X2)  -  e  '  g^.  Also,  by  (2.13), 

IhMxj)!  =  /hh"  +  0(e)  =  o(e’^^)  . 


Therefore  by  (2.18) 

0  <  -G'(X2)  •=  0(e’^^)  . 

By  Lemma  2.1  the  bound  on  H'*'  depends  only  on  the  bound  for  G/H  which  is  independent 
of  M.  Therefore  H'*'  can  be  estimated  independently  of  u.  This  proves  the  lemma. 

Now  we  can  proceed  as  in  Section  3.  For  x  >  Xj  we  introduce  new  variables 
h(x)  =  H(x)/e^'^^,  g(x)  »  G(x)/e^^^,  x  -  X2  =  e'^^x,  x  >  0  . 


Then  we  obtain  the  equations  (3.8) 


(4.7) 


h  +  hh  +  -^  -  h^) 

g  +  hg  -  gh  =  0 


with  boundary  conditions 
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(4.U) 


h(0)  «  hy,  h(0)  »  ,  h(0)  «  hj,  g(0)  ■  gp,  g(0)  •  g, 


Also 

(4.9) 

’h*  (0)  = 

'■'(Xj)  «  0(e*^^),  h'*'(x2)  is  bounded  independently  of  u 

and 

(4.10) 

5o/''o  ■ 

5l  " 

T|  =  G(x, )/H(x^)  +  0(1/h^/^)  «  T  +  o(e  t  hp^/^) 

T^h^  +  0(1/hp)  -  Thj  +  0(e  +  l/hp)  . 

(4.9)  ana  (4.7)  snow  tliat 

(4.11)  2hQh2  -  -  0(e’''^),  i.e.  h,  =  -  +  0(e’‘^^)  , 

ana  the  uniform  boundedness  of  ri*''(x)  for  4  x  4  X2  sh^ws  t)uit 

(4.12)  lim  h2  ~  H2  * 

e*0 

We  can  solve  Che  reduced  problem 

•  i  •  ••  i  •• 

(4.13)  h  +  hh  -  ^  «  0,  h(0)  »  hp,  h(0)  -  h, ,  h(0)  -  h2 

•«  •taai 

(4.14)  g  +  hg  -  gh  *  0,  g(0)  =  g(0)  •  g^  . 

By  (4.11)  the  solution  of  (4.13)  is  in  any  finite  interval  0  4  x  4  x^  of  the  form 

(4.15)  h  = -j  hjCic  -  Xj)^  +  0(e’''^) 
where  x^  is  determined  by 


’  '2 
2  ^*3 


+  o<c’/^). 


/5. 


“  1  2 

Replacing  h  in  (4.14)  by  —  h  'x  -  x,)  and  introducing  a  new  variable 

V3  ^  ^ 

t  -  (^  ‘'2^  ■  *3’ 


gives  us 

(4.16)  d^g/dt^  +  t^dg/dC  -  2Cg  »  0,  ^  =  -[-j  h^)  *3  "  * 


L  (1 


1/3 


1/2 


*o(c’/^ 


with  boundary  conditions 


(4.17) 


g(ep) 


dg(Cp)/d5 


1  -V3 

^1 


By  (4.10)  and  (4.11) 
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(4.1«) 


^0  ,1  -3/2 


1 

a5(«o)/d«  (2  *>2)  g,  ,/3 

+  0(e  '  ) 


2«„ 


-./I^ 


,  1/3 

,  0(eV^ 

0  2 


1/3 


(■J  1*2^  ^  *  0<h"3/2  ^  e1/3)  _ 


Thus  by  the  Appendix 


J/3 


5(t)  =  (-Jhj)  Tg^(e)  +  0(h’^''^  + 
where  g^(0  Is  monotone  decaying  with 


q,(C)  9,(5) 

lim - -  =1,  9,(0)  <  0,  lim  "■  - 

6*--  r  5*+-  5^ 


.|t|  -  T  . 


For  the  original  equations  (4.7)  a  standard  perturbation  analysis  gives  us  in  any  tinite 
interval  0  <  x  4  x  , 


(4.19) 


li(x)  -  j  h2(x  -  x^)^  +  0(e'^^)  , 


1/3 


g(x)  =  (j  h^)  Tg^  ([y  h^  J’'^^(x  -  x^))  +  0(h"^^^  +  e’''^)  . 

Thus  we  can  shoot  through  the  turning  point.  In  particular,  we  can  choose  x^  so  large 

that  the  conditions  of  Theorem  3.1  are  satisfied.  Furthermore 

g(x  )/h(x  )  =  TT  +  u(h-^^^  +  c^^^)  . 

4  4  0 

Thus  we  can  continue  the  solution  of  our  problem  to  the  next  turning  point 

X  -  2»(-^  +  — — )  where  we  can  repeat  the  process. 

In  the  next  section  we  need 

Lemma  4.2:  For  e  sufficiently  small  H"'(x)  can  be  estimated  independently  of  u. 
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Prooti  Away  tcon  tha  turning  pointa  this  la  clear  because  by  Lemna  2. l  a  bound  for 


H"'(x)  depends  only  on  a  bound  for  G/H  ^  t.  To  estimate  H' "  for  0  <  x  <  we  can 
2/3  2 

n«glect  the  tern  €  '  g  in  (4«7)«  Differentiating  gives  us 

X 

-  /  h(«)de 

•  •••  •••  •••v  r>  ••• 

h  +  h  h  ■  0,  h  (x)  >•  e  h  (0)  . 

*  *  *  t  /3  - 

We  know  that  h  (0)/e  ^  is  bounded  independently  of  m  and  therefore  the  same  is  true 
•  • *  ^  1/3 

for  h  (x)  because  h(C)  >  •‘0{e  '  )«  This  proves  the  lemma* 


5.  Existence  of  Solutions 


Let  n  >  1,gQ,g^  be  given  with 


sign  g,  =  (-1 )" 


where  T  is  the  constant  given  in  (A.6).  Let 

(!>.2a)  h^  =  |(g,/l100Tg^(0)l)^l 

where  g^(x)  is  the  function  described  by  (A.1)  with  6  s  i  and  (A. 2a).  Let 
(5.2b)  hj  -  |(200g^/ttg^(0)l)^|  +  1  . 

Let  Q  be  the  bound  on  |h'''(x)|  deternined  in  Lemma  4.2.  Let 
(5.3)  8  «  Q  +  1  . 

For  the  remainder  of  this  section  we  require  that 


|ti|  <  6,  i.e.  |v|  <  Be,  =  1  , 


(5.4b)  iia  ‘  >'2  *  ^2  • 

Let  (G(x,e)ii,h2)>H(x,eiU,h2))  the  solution  of  (2.1),  (2.2)  which  arise  from  the 

pair  <g(4,e;u,h2),h(4,EtM,h2))  which  satisfy  the  boundary  value  problem  (3.8),  (3.9). 

The  constant  a  is  fixed  with  a  >>  1. 

By  the  arguments  developed  in  the  previous  sections  we  obtain  the  following  results. 
Theorem  5.1;  Let  n  >  1  be  a  fixed  integer.  If  a  is  chosen  large  enough  and  e  small 
then  (  G(x,C)  |i,h^ )  ,H(x,e)  u.h^))  exist  on  an  interval  (0,Bl  whose  length  B  is  of  order 
1.  The  function  G(x,e;|i,h2)  has  at  least  n  zeros 
(5.5a)  0  <  X,  <  Xj  <  •••  Xj^  <  B  . 

The  function  IKxjejU.h^)  has  at  least  n  relative  minima 
(5.5b)  0  <  y,  <  y2  <  •••  <  y^j  <  B  . 

It  gjj  <  0  then  G(x,e;|i,h2)  has  another  zero  Xq  with  0  <  Xq  <  x^.  These  numbers 


satisfy 


0  <  Xj  -  y^  4  0(1)  as  e  ♦  0,  j  =  1,2,< 
X,  =  2»  +  o( 1 )  , 

X..,  -  X,  =  (2x)/|t1^‘’  +  0(1),  j  =  1,2, 


Also 


(S.bd)  H"(x  ,e,u,h_)  »  h.,  +  o(1)  . 

n  z  ^ 

On  the  interior  of  the  interval  j>1,2,...,n  we  have 

h  ^  . 

(5.7a)  HCx.ejii.h^)  »  ^^^(j-T)  ”  cos|t1^  (x  -  x^_,))  +  o(1) 

while 

(5.7b)  ^(x.eiM.h^)  «  (T)^"’H(x,e;p,h2)  +  oO)  . 

It  u  a  lie  then 

(5.8a)  H(x  ,e>8,h  )  >  0  , 

n  z 

and,  if  ii  «  -$e  then 

(5.8b)  H(x  ,e;-$,h,)  <  0  . 

n  2 

It  h2  »  ^2*  then 

(5.9a)  G(x^,e,M,h2)  ”  ^2^  I ""'g,  ( 0)  ]e^^^  +  o(e^'^^)  >  . 

similarly,  if  h2  ■*  2^2' 

(5.9b)  G(x^,e,w,h2)  •  [(j  hj)  |T|"”’g^  (0)  +  o(e*'^^)  <  9,6^^^  . 

Corolla ryi  There  is  a  choice  of  v.h^  which  satisfy  (5.4a),  (5.4b)  and 

(5.10a)  H(x  ,e>u,h_)  •  0 

n  z 

•  2/3 

G(x  , E;u,h_)  »  g.e  . 
n  2  1 

Proof!  Let  e  >  0  be  so  small  that  Theorem  5.1  holds.  Consider  the  mapping 

(5.11a)  (M,h  )  ♦  <G(x  ,e)U,h  ),H(x  ,efii,h,)>  . 

2  n  2  n  2 

With 

(5.11b)  |p|  <8,  ilj  *  »*2  ‘  • 

Since  H"(x^,e,U!h2)  "  hj  >  0  the  implicit  function  theorem  shows  that  x^^  defined  by 
the  n'th  mint 

HMx  ,e;li,h,)  -  0 
n  2 

is  d  continuuus  function  ot  Thus,  the  mappinq  (5s11a),  (b.llb)  is  a  continuous 

niappi  ng. 
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The  properties  (5.5a),  (5.5b),  (5.(>a),  (5.bb)  together  with  an  elementary  "degree" 
theory  argument  see  (2UI  shows  that  there  is  a  solution  of  (5.7a),  (5.7b).  To  see  this  we 
consider  the  homotropy 


*  3  ^2  *  ^2  2y3 

^tCx^.ElU.h^)  -  tgja  - ^)]e  '  +  (1  -  t)G(x^,e;U,h2),  0  <  t  <  1  , 

h_  -  h^ 

2  —2 


H^(x  ,e;U,h..)  =  tn  +  (1  -  t)H(x  ,ejU,h.),  0  <  t  <  1  . 
t  n  2  n  2 

As  t  varies  from  0  to  1  the  inequalities  (5.5a),  (5.5b),  (5.6a),  (5.6b)  continue  to 
hold.  Thus,  throughout  the  homotropy,  there  is  no  solution  on  the  boundary  of  the  region 
described  by  (5.11b).  For  t  =  1  the  equations  read 


(5.12a)  [2  -  I  (^^2 - 2)]g  .  , 

*  ^  *  i2.2 

(5.12b)  U  »  0  . 

There  is  a  unique  solution, 

**2  “  T  ^2  *  1  —2'  >*  “  0  • 

Thus,  there  is  a  solution  h2(t),M(t)  for  every  t  6  [0,11.  In  particular,  there  is  a 
solution  for  t  ..  1  and  our  problem  has  a  solution  (see  [20]). 

This  corollary  implies  the  truth  of  Theorem  I.  Theorem  II  and  Theorem  III  follow  as 
indicated  in  the  Introduction. 
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Appendix 


In  thl«  Appendix  we  are  concerned  with  the  equation 
(A.1)  q"  +  Sx%*  -  2Sxq  «  0,  fi  >  0  • 

In  factf  we  need  only  consider  the  case  6-1.  For,  If  q(x;1)  Is  a  solution  of 
(A.1)  -  with  6  -  1  -  then,  for  any  6  >  0,  a  direct  calculation  show  that  the  function 

Y(x>«)  -  q(«’^^x,1) 

la  a  solution  of  (A.1)  with  this  value  of  6. 

Our  first  concern  la  with  the  asyeptotlc  behavior  of  solutions  q(C)(-  q(€,1))  as 

C  *  ±". 

A  alsple  calculation  uslnq  the  Uouvl  lle-Sreen  (or  HXBJ)  approximation  (see  chapter  6 
of  (14))  leads  to  the  followlnq  results. 

Case  1 1  As  X  there  are  two  linearly  Independent  solutions  q^(x),a2(x) 

(A.2a)  q,(x)  ~  x^,  X  ♦  -«  , 

3 

(A.3b)  q2(*)  "*  *”  exp[”^]«  X  ♦  . 

nitts,  there  is  a  unique  function,  q^(C),  <d)lch  satisfies  (A.1)  with  6  -  1  and 
(A.3)  q,(x)/x*  ♦  1,  X  . 

Furthermore,  a  more  careful  asymptotic  expansion  of  g^(x),  e.q. ,  using  the  methods 
described  In  (25,  pp.  52-61)  yields, 

(A.4b)  g,(x)  ~  x^(1  2/3x^),  X  ♦  , 

(A.4b)  g^(x)  ~  2X(1  +  2/3x*)(1  -  1/x^),  X  *  , 

and 

(A. 4c)  q^(x)  ~  2,  X  ♦  . 

Case  2«  Ihe  same  calculations  show  that*  as  x  *  there  are  two  linearly  independent 
functions  ^^(x),^2(a)  which  satisfy  (A.1),  with  6-1,  and 

3 

(A. 5a)  y>,(x)  ~  x“>xp[- ^],  x  ♦  ■*»  , 

(A.5h)  '/>2(x)  X*,  X  ♦  . 

Since  the  function  q^(x)  which  Is  characterized  by  (A.3)  can  be  written  as  a  linear 


I 


combination  of  ^^(x)  and  that  there  is  a  unique  constant,  call  it  T, 


such  that 


Lim  g, (x)/x  =  i  . 

x-»+» 


Lemma  A.  1 ;  Let  i^^Cx)  be  the  solution  of  (A.1)  -  with  6  -  1  -  described  by  (A. 5a).  Then 


(0)  *  0  . 


Proof:  Suppose  (A. 7)  is  false.  Then 


A  «  \^i’(0)  *  0  . 

Suppose  A  >  0.  An  easy  argument  based  on  the  maximum  principle  -  or  based  on  the 


identity 


shows  that 


The  identity 


and  the  fact 


implies  that 


('P*exp[^j)}  =  2x.p^exp[^] 


i(>,’(x)  >  0,  0  <  X  <  -  . 

^,(x)  >0,  0  <  X  <  "  . 


3  3 

—  ('#>;|exp[^]}  »  2'^^exp['2^] 


^j(O)  =  0 


.p;[(x)  >0,  0  <  X  < 


(x)  >  Ax 


which  contradicts  (A. 5a).  If  A  <  0  we  apply  the  above  argument  to  -<^^(x). 

We  now  turn  to  a  more  detailed  discussion  of  the  function  q.|(x). 

Theorem  A;  Let  gj(x)  be  the  function  which  satisfies  (A.1)  -  with  5=1-  and  (A. 3). 


(A. lUa) 


(A. 10b) 


q,'(x)  <  0, 
xq?  <  0. 


Let  g  uenote  the  unique  point  at  which  q^(x)  vanishes,  i.e.. 
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(A. 11a)  q^Cq)  »  0  • 

Then 

(A. 11b)  -<2/3)’/^  <  q  <  0  . 

Finally,  the  number  T  of  (A«6)  satisfies 

(A. 12)  T  <  0  . 

Proof;  For  neqative  values  x  «  -  R  <  0  we  have 

q^(x)  >  0,  qj(x)  <0,  X  <<  -R  <  0  . 

Since  a^(x)  satisfies  (A.1),  q^(x)  cannot  have  a  positive  relative  minimum  on  the 

Interval  (-'•,0).  Thus,  either 

(A. 13)  q,(x)  >0,  -»  <  X  <  0  , 

or  there  is  a  first  point  q  <  0  at  which  (A. 11a)  holds.  Suppose  (A. 13)  holds.  ‘Aien 
(A. 9)  shows  that 

q^(x)  >  0,  -»  <  X  <  0  . 

Let  X.,  <<  -1  and  x.,  <  x  <  0.  Tnen 


(A. 14) 


0  <  q,(x)/x^  ■  q,fx,)/xj  -  /  (q^(t)l/t^dt  <  q,(x,)/x*  . 


Let  x.|  ♦  Then 


0  <  q^(x)/x^  <  lim  q,{x,)/x. 


1 


1  . 


Thus,  if  (A. 13)  holds. 


But  then 


and 


0  <  q,{x)  <  x*,  -•»  <  X  <  0  . 

q,  (0)  »  q^(0)  «  0 


q,(x)  so. 

Since  this  is  impossible,  there  is  a  first  point  q  <  0  at  which  (A. 11a)  holds.  Moreover 
(A. 15a)  q{(x)  <0,  -••<x<q<0, 

(A.ISb)  q"(x)  >0,  -“rxtqrO. 
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we  now  estimate  y  trom  below*  We  have 

<jj"(g)  »  2q,(q>  -  q^gl[(g)  <  0  . 

Let  X,  <  q  be  any  point  such  that 

(A. lb)  q'"(x)  <  0,  X,  <  X  <  q  . 

Then 

0  <  qi(g)  <  g^Cx)  <  g^(x, ),  x,  <  x  <  g  <  O  . 
Returning  to  (A. 14)  we  have 


g,(g)/g^  =  O  =  q,(x^)/x^  -  /  [g"(t)/t^)dt  . 


That  is 


q  q 

U  <  q,(x,)/x^  =  /  [q^(t)/t^ldt  <  q',’(x,)  / 


[1/t*]dt  . 


Thus 

(A.17)  0  <  g,ix,)/x^  <  [g!|(x,)/3)((1/|g|^)  -  (1/|x,p))  <  g;[(x,  )/(3  |g|  ^)  . 

Let  X,  ♦  under  the  condition  that  (A. 1b)  hold.  Either  x,  ♦  x,  a  finite  point  at  which 

qJ"<X,)  =  0 

or  X,  ♦  In  either  case 

2g, (X, )/xJ  ♦  qj (X, )  . 

Thus,  (A.17)  yields 

0  <  |g|^  <  2/3  . 

Thus  we  have  proven  that  (A. 11b)  holds  for  the  "first"  zero  of  g^(x). 

Our  next  tas)(  is  to  extend  the  range  of  the  inequalities  (A.iba),  (A. 15b)  to  the 
larger  interval  (-".O). 


Let 

6  » 

1  and 

let 

7,  (x),Y2(x) 

be  the 

special  solutions  of  (A 

. 1 )  which 

also 

(A.iba) 

Y,(0) 

=  0, 

y;(0)  =  -1  , 

(A. lUb) 

Y^CO) 

-  -1, 

YJ(0)  -  0  . 

Let  tj 

and 

r^  be 

the 

smallest  (in 

absolute  value)  negative  zeros 

of  Y,  (x) 

and 

Y^(x)  rf'S|)ectivel  V" 

It 

IS  an  easy  matter  to 

obtain  intinite  series 

soliitions 

tor 
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y^(x),y2(x)  and  sea  that 

[rj  >  1.  |r2l  >  1  . 

These  estimates,  together  with  (A. lib)  and  the  oscillation  theorems  (see  (16,  page  421) 
show  that 

g,(x)  <0,  g  <  X  <  0  . 

Thus 

(A. 19)  g,(0)  <  0  . 

l,et 


g,(x)  »  d,.',(x)  ♦  d2y2(x)  . 

Since  g^(g)  0  and  y^(g)  >  0,  y2(g)  <  0  we  see  that  d^  and  d2  are  ot  the  same 

sign.  But  (A, 19)  gives 


Thus, 


d2(-l>  -  9,(0)  <  0  . 


and 


d,  >  0,  dj  >  0  , 


g,*(0)  -  <  0  . 

But,  (A.1)  implies  that  g, (x)  cannot  have  a  negative  relative  maximum  in  the  interval 
(g,0).  Since  g{(g)  <  0  and  g|(0)  <  0  we  see  that 

g{(x)  <0,  g  4  X  4  0  . 

On  this  interval  we  have 


g^  »  2xg,  -  x^g^  >0,  g  4  x  <  0  . 

Hence,  »re  have  the  inequalities  (A.IOa)  and  (A.  10b)  on  the  interval  (-",0].  The 
completion  of  the  proof  now  follows  from  the  initial  conditions  g,(0),gf(0),  the  maximum 
principle  and  the  identity  (A,9). 

Remark;  Since  the  theorem  holds  the  quantities  g  and  x  can  be  determined  -  to  any 
desired  accuracy  -  by  numerical  computations.  Results  of  Jerry  Browning  of  NCAR  indicate 
that 


t  '  -2  , 


g  -.91  . 
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